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Abstract 

The allocation problem for a Poisson point process is to find a way 
to partition the space to parts of equal size, and assign the parts to 
the configuration points in a measurable, "deterministic" (equivariant) 
way. The goal is to make the diameter of the part assigned to a 
configuration point have fast decay. We present an algorithm for d > 
3, that achieves an 0(exp(— cR d )) tail, which is optimal up to c. This 
improves the best previously known allocation rule, the gravitational 
allocation, which has exp(— cR 1+ °^) tail. 

1 Introduction 

We prove the existence of an allocation rule of the following (optimal) tail. 

Theorem 1.1. For d > 3 there exists a c > and an allocation rule i/j for 
the Poisson point process such that 

P[diam({0} U V(0)) > R | E u] < exp(-cR d ). 

From now on, we always assume that d > 3. Consider Poisson point 
process of intensity 1 in M. d . An allocation rule is, informally, a rule to 
partition the space into finite pieces, each of equal measure, and to give a 
bijection between the Poisson points and the pieces, in such a way that the 
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partition is (deterministically) determined only by the point configuration, 
and the rule is equivariant. That is, for any configuration point £, for any 
e > there is an R(e) such that with probability at least 1 — e a cell of £ can be 
determined by looking at the configuration points in the i?(e)-neighborhood 
of £, and this cell has Hausdorff distance at most e from the cell assigned to £ 
by the allocation rule. This latter property is the requirement of measurabilty 
of the allocation rule. Even more informally, in the allocation problem we 
want to divide land to a set of farmers scattered in space, and in such a way 
that each farmer can determine his own land up to a small error, by looking 
around a big enough neighborhood, and this rule is the same for everybody. 
One wants to find an allocation rule where the maximal distance a farmer 
has to walk (fly) to arrive to any point of his land is minimal. That is, we 
want the tail 

P[diam({0} U ^(0)) > R | e u] 

decay as fast in R as possible, where ^(0) denotes the cell of (conditioned 
on having a configuration point in 0). 

The allocation problem was first studied in a finite setup, where finitely 
many points are distributed uniformly and independently in a box. Here, 
of course, the requirement of equivariance is meaningless. We will present 
later a variant of the algorithm by Ajtai, Komlos and Tusnady p], which 
was the starting point of all later methods. For n uniformly independently 
distributed points in a cube of volume n, it was proved ([I]) that the average 
diameter of an allocation cell is log 1//2 n for d = 2 and finite for d > 3; 
and precise rates were determined subsequently (see [9] for details and the 
sharpest results). Interest in the infinite setup originated from the existence 
of a shift-coupling between a point process and its Palm version, see e.g. 
[5], [B]. It was first asked in [5], what the optimal tail of a nonrandomized 
allocation rule is. In the same paper it was shown that constructing extra 
head rules and invariant allocations are equivalent, furthermore a randomized 
invariant allocation rule of optimal tail decay was presented. Deterministic 
constructions have been an active research area ever since, allocation rules 
with respect to additional conditions have been subject of recent analysis 
(e.g. stability 0], connectedness [8]). The best previously known allocation 
rule for d > 3 was the gravitational allocation , investigated in [2], [3], where 
the tail is P[X > R] = exp(—R 1+ °^). The assumption d > 3 is necessary 
for an exponential tail by [6J, where E[X d / 2 ] = oo is proved for d = 1,2. 
The brief reason for the drastic change of behaviour from dimension 3 is that 
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here the isoperimetric function of M. d becomes larger in magnitude than the 
deviation of the number of points in a ball. For the existence of optimal 
allocations with respect to other quantities (e.g., the integral of the distance 
of a center from the points of the cell) and connections to optimal transport, 
see [TJ. 

Definition (Allocation) Let Q be the set of discrete sets of points in 
M. d . The points of an u G Q are called centers. An allocation is a function 
ip u = i/j : u -> L 1 (R d ) such that 

(1) For every £ G w, t/>(£) takes values from {0, 1}. 

(2) For Lebesgue almost every x G M d there is exactly one (€w such that 

= 1. 

(3) For every £ G w, / Rd ^(f)dac = L 

Without assumption (3) we call ip a wea/s allocation. We say that ^ is a 
weak fractional allocation, if the following (1') and (2') hold, and a, fractional 
allocation if (1'), (2') and (3) hold: 

(1') For every £ G u, ■?/>(£) takes values from [0, 1]. 

(2') For Lebesgue almost every x G M. d , Y^tew ip(£,)( x ) = 1- 

Sometimes it will be natural to think about a (fractional) allocation ip : 
uj — > L 1 {R d ) as a family {"0(0 : ^ G w} of functions. 

Consider an allocation, weak allocation or fractional allocation ip, and 
£ G w. By a slight abuse of notation we define the diameter of as 
max{|a; — y\ : x,y E supp(^(£)) U {£}}. (Note that this is slightly different 
from the usual diameter of the set t/>(£)). We denote it as diam(^(£))- We will 
also use diameter in the usual context. There will be no danger of ambiguity; 
what matters is that when we are talking about the diameter of a cell, we 
always include the corresponding center to the cell, for simplicity of notation. 

Definition (Allocation rule) Let V be the Poisson point process of 
intensity 1 in IR d . An allocation rule is a mapping u i— > that is defined for 
P-almost every uj G Q, measurable (with respect to the relevant a-algebras), 
and such that: (i) almost surely t/> u is an allocation, and (ii) the mapping 
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co i — y ip^j is translation-equivariant, i.e. for any w 6 fi, ( 6 w and x,y E M. d 
we have 

+ + x) = V«(f)(2/)- 
Define u>ea£;, fractional and u>eafc fractional allocation rules analogously. 

We want to define an allocation rule in such a way that the probability 
P[diam( , 0(O)) > R | G u] decays as fast as possible. By translation invari- 
ance we could have taken any other fixed £ G M. d instead of 0, and the tail 
would be the same. We will use notation u — £ for the translate of u by — £ 
(so (£w implies G w — £). 

Define Q' := {w U {0} : u £ Q}, and P' as the Palm version of the 
probability measure of the Poisson point process (and E' the corresponding 
expectation). It is well known that P'[{w U {0} : u G E}] = P[E] for every 
measurable E inside Q. 

In several papers ([2], [I]) a different problem is in the focus of interest. 
Suppose we have an allocation rule uj ipu, this defines a unique config- 
uration point £ — £o( w ) with G V'w(^o) almost surely. The objective in 
this setup is again to obtain a rule with optimal tail bound, this time for the 
random variable |£o|- Our setup is stronger than this one and applies to this 
optimization problem, since 

P(|£„| > R) < P(max xe ^ (Co) |s - £ Q \ > R) < P(diam(^(£ )) > R) = 

= P(diam(^_ ?o (0)) >R) = P'(diam(Vv(0)) > R), (1.1) 

where the last equation follows from Theorem 13 in [6], which verifies that 
£o is a nonrandomized extra head scheme. In fact, imagine an allocation rule 
when most cells consist of ball-like pieces of almost unit volume, and a small 
extra piece far away from this one. For such an allocation, the quantity on 
the left of fll.lD decays fast, while on the right we may have slow decay. 

Let us sketch the (surprisingly simple) construction briefly, before going 
into the details. First, we will define a weak allocation scheme associated 
to a i) £ M. d , which will not necessarily be invariant. This will be based 
on a straightforward generalization of the algorithm of Ajtai, Komlos and 
Tusnady (AKT algorithm), [1], which assigns equal volume to each of n 
points in a box B of volume n, in such a way that the pieces partition B. 
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Furthermore, if the points are scattered uniformly and independently (which 
is the same as the Poisson case, conditioned on that there are n points in B), 
then the average diameter of a piece has asymptotically same tail behavior as 
in Theorem 1 (regardless of n). The algorithm proceeds as follows. Subdivide 
B to dyadic cubes, until each cube contains at most one configuration point. 
Then, starting from the level of smallest cubes, match cubes with each other 
in the simplest possible way (so that the resulting pairs are all translates of 
one another), and translate the "wall" between each pair in an affine way, 
together with the configuration points inside, and in such a way that the 
volumes after the translation are proportional to the number of points in 
the respective parts. (See Figure [TJ and Section 2 for more details.) Call 
the procedure again for the next level (with cuboids of edge lengths c and 
2c), and repeat, until reaching the biggest cube. The resulting cells all have 
equal volumes, and we assign each such cell to the original configuration 
point whose image after the repeated translations is in the cell. 

The AKT weak allocation in an infinite setup will work as follows. (This 
will not be a weak allocation rule yet, because equivariance fails.) Assume 
for a second that the partition v + 7L d is such that each class (cube) contains 
at most one configuration point from u. Then we can repeatedly apply the 
previous procedure, with cubes in the i'th stage coming from the partition 
v + 2 l 1* d , as % = 1, 2, . . . , n. Say that we ran the AKT(v) algorithm up to 
stage n for the result of this weak allocation. We can remove the assumption 
that each class of v + Z d contains at most 1 point from oj: simply subdivide 
diadically those unit cubes B where \B fl u\ > 2 repeatedly, and for each B 
start the AKT procedure from the lowest level of this subdivision (but still 
run it up to v + 2 n Z d ). 

Next we will define an equivariant weak fractional allocation from the 
above allocations. Fix u' G Q', and define the following deterministic func- 
tions of oj' (we will often hide the dependence on u' for simplicity). For 
v G M. d , let f£ n = f Vt n be the indicator function for the cell of the origin 
constructed by AKT(t>) run up to stage n. Note that 



since the sequence of cubes up to stage n given by such a u and v is the same. 
Now, define 



v, n Ju,n 



whenever u — v G 2 n Z d 



(1.2) 




(1.3) 
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a function from M. d to [0,1]. This is well defined (the integral exists), by 
Lemma [3.21 below. We will prove that the f n have an L 1 limit /, and that 
it is a function of integral 1. From the properties of the AKT construction 
(Lemma 12.61 below) we will get that the diameter of this limit has the tail 
that we want. We will conclude that the map ip : £ H> f u ~* (f G u G fi) 
defines a fractional allocation rule with the desired tail. 

Finally, we will define an allocation from the above fractional allocation. 
It will be such that the support of the allocation function corresponding to a 
configuration point £ is contained in the support of the fractional allocation 
function f^^^ corresponding to the point. This step (Lemma 13.71) will be 
based on the fact that almost every point of M. d is contained in the support 
of only finitely many functions of the fractional allocation defined above. 

In the next section we give a summary of the AKT construction, pre- 
senting the generalized version that we are using, together with the neces- 
sary bounds for the expected diameter and concentration around it. The 
third section continues with the sequence of invariant weak allocations, the 
limiting fractional allocation, and finally, the allocation rule that satisfies 
Theorem 11.11 

2 The generalized AKT algorithm, bounds 

Our main tool is the so-called Ajtai-Komlos-Tusnady transport scheme which 
leads to allocation in a fixed cube. We will generalize the approach for the 
infinite setting. Take the Poisson point process and the fixed cube [0,2^]^ 
for some N G N (see Figure [1]). Assume that there is at most 1 configuration 
point in each dyadic cube in {(ki, . . . , k d ) + [0, l] d : < ki < 2 N — 1 for 1 < 
i < d}. 

For the first step, consider the cuboids {(ki, . . . , k d ) + [0, l] d_1 x [0, 2] : < 
h < 2 N - 1 for 1 < i < d - 1 and < k d < 2 N ^ - 1} and take the bisector 
hyperplane orthogonal to the GPth axis of each of these cuboids, and move 
them along the d'th axis until the volume ratio of the two halves gets equal 
to the ratio of configuration points inside them, and transform the respective 
interiors affinely (do this separately for each of the cuboids). Now proceed 
to the second step, consider the cuboids {(ki, . . . ,k d ) + [0, l] d ~ 2 x [0,2] 2 : 
< ki < 2 N - 1 for 1 < i < d - 1 and < k d ^, k d < 2 N ' X - 1} and take 
their bisector hyperplanes orthogonal to the (d — l)'th axis, and move them 
like before, for each cuboid separately, interiors transform affinely Continue 
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Figure 1.1: The Ajtai-Komlos-Tusnady transport scheme in the interior of 
a fixed cube of sidelength 4. The translates of the centers indicate how the 
cells will be assigned to the original centers. 



this procedure of transformation of cuboids always of double size until the 
cube [0, 2^]^ is reached (at this last step we have only one cube to consider), 
the result is a partition of [0, 2 N ] d , where the cells are transformations of the 
initial cubes with exactly one translated configuration point inside them, and 
all cells have equal volume. The cells can be assigned to the untranslated 
configuration points through the transformed ones inside them. Let us call 
each transformation a step, as above, and the (i-tuple of the first d steps, the 
c?-tuple of the second d steps etc. stages. 

Some care has to be taken when our assumption does not hold, i.e. there 
are initial cubes with more than one configuration point, call them "bad". In 
that case consider each "bad" initial cube separately. Subdivide them one- 
by-one dyadically until all the subcubes of sidelength 2~ k contain at most one 
configuration point. This k exists almost surely for each unit cube (but we 
will hide dependence on the unit cube), so the following is well-defined. First, 
run the natural scaled version of the previous algorithm in the "bad" cubes 
of sidelength 1. This algorithm will equipartition these cubes. Then run 
the previous algorithm as before, for each cube of sidelength 1 (transforming 
cells within "bad" cells affinely) with considering all walls inside the bad 
cubes. This will result, similarly as before, an equipartition of [0,2^]^ into 
cells assigned to the original configuration points. For notational purposes 
index the initial steps in the "bad" cubes with —kd+ 1, — kd + 2, . . . ,0. 

Given a v G M. d , we can run the above transport scheme for each 2^ x 
. . . x 2 box of v + 2^!^. The resulting partition of M. d is what we called 
the AKT(v) run up to stage N. (Of course, it is by no means an invariant 
partition, since it depends on v.) We define the v-partition sequence as 
the series of coarser and coarser partitions determined by v + Z d ,v + 2Z d , 
...,v + 2 ir L d , . . . with v e R d . 

In each stage there is exactly 1 step along each of the d axis. Also, steps 
along different axis are independent. Therefore if we want obtain an upper 
bound on the total movement of a point x, steps along different axis can be 
treated separately. 

Condition now on 6 w, and let C v>n be the cell assigned to by AKT(u) 
run up to stage n. Let f v>n be the indicator function of C v>n . According to the 
sketch of our construction, in order to prove Theorem we have to study 
first the tail behavior of diam(C„ jn U{0}). After n stages, we can think of the 
cell C v>n as the transformed initial cube of after nd afline transformations 
((n + k)d respectively, if was in a bad cube). Hence, we can bound the 
diameter of C v>n U {0} by bounding first the total shift of after n stages, 
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and bounding the length of the edges of C V)H . 

The next lemma is standard, we prove it only for completeness. 

Lemma 2.1. Let X be a random variable with Poisson distribution of mean 
A. Then ifO < p < 2 then 

P(|X- A| > Ap) < 2exp( — 

Proof. Note that the moment generating function of X is: 
M(t) = E(exp(tX)) = exp(A(e* - 1)) te R. 

For one side 

P(X - A > Ap) = P(exp(tX) > exp(At(p + 1))) < 

E(exp(tX)) exp(A(e* - 1)) 
< exp(At(p+ 1)) ~ exp(At(p+ 1))' 

where t > 0, and we used Markov's inequality. Now, for < t < 1 we have 
e* < 1 + t + t 2 , so 

P(A ' " A > A "» < M + l)) = exp(A((2 + «) " Ai( " 

The last expression is minimized by t = p/2, so 

P(X-A>Ap)<exp(-^). 

For the other bound: 

P(A - X > Ap) = P(exp(-tX) > exp(-At(l - p))) < 

E(exp(-Xt)) exp(A(e- t - 1)) 
< exp(-At(l - p)) ~ exp(-At(l -p))' 
where t > 0, and we used again Markov's inequality. Now, for < t we have 
e - * < 1 - t + t 2 /2, so 

P(A - X > Ap) < ^^"^ = exp(A(t 2 /2 -t) + Xt(l - p)). 

The last expression is minimized by t = p, so 

P(A-X > Ap) <exp(-^— ). 
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□ 

For a measurable subset B C M. d , let N(B) denote the number of config- 
uration points of the Poisson point process in B. The next lemma delivers 
the desired information about the displacement of configuration points. 

Lemma 2.2. There exist c,C > such that for all R > there exists an 
event Er such that 

(i) 1 - P'(E R ) < Cexp(-cR d ). 

(ii) On Er, for every n, the total movement ofO in the AKT(v) run up to 
stage n is at most R, for all v G M. d . 

(Hi) On Er, the total length of the shifts of in the i'th stage is at most 
c\ = c'^R), and {c'A- is absolutely summable. 

Proof. As we do not wish to optimize on the constants, we may prove the 
lemma only for R > R for some suitable Rq > 0. That implies the lemma 
for every R, possibly with different constants. 

Fix Del d and R > Ro, Ro will be determined later. We remark that Er 
will depend on this arbitrarily chosen reference point v , therefore it should be 
the same for all R. Our analysis will loosely follow the argument of Talagrand 
and Yukich j9], although we are allowed to use less sophisticated methods (i.e. 
triangle inequality), as sizes of displacements induced by points of increasing 
distance decay much more rapidly in d > 3 than in d = 2. Instead, we have 
to achieve uniform bounds regarding the f-partition sequences. 

As, informally speaking, shifts of in each direction at one particular 
stage depend only on the number of configuration points in cuboids, and 
these cuboids have comparable volumes in the particular stage, it is enough to 
consider only shifts along the first axis, the others can be treated completely 
analogously with essentially the same bounds on the length of displacement 
of points, with the same probabilities. For bounding the size of total shift it 
is then enough to use triangle-inequality. 

We will define an event Ar = Ar^, that delivers bounds on the number 
of configuration points in both halves of finitely many cubes for each n > r 
(ro to be determinded later) ("left" and "right") with the following three 
properties: the cube has corners in G n (v) := v + 2~ n Z d , sidelength 2 n ± 2~" 
and contains the origin. The second parameter of the event Ar^ refers to 
the direction of the shifts in consideration. 
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This conditioning will imply bounds on the number of configuration points 
in all cubes of sidelength 2 n containing (i.e. not necessarily only for those 
with vertices in G n {v)). 

Let r be such that j < 2 r ° < -|; and for n = r + 1, r + 2, . . . , mi = 
1, . . . , (2 2n - l) d and m 2 = 1, . . . , (2 2n + l) d let Q n>mi denote the mi'th cube 
with vertices in G n (v) of sidelength 2 n — 2~ n containing the origin (we assume 
that the cubes are numbered in some arbitrary way), U njTni its left side, V n>mi 

its right side; and Q n ,m 2 denote the m 2 'th cube with vertices in G n (v) of 
sidelength 2 n + 2~ n containing the origin, U n>rn2 and V n ^ m2 denote the two 
halves respectively. Let Ar be the event, when all the following hold: 

/on o — n\d /on o — n\d 

AT 1 ■■ \N(Un, mi \ {0}) - 1 2 } I < 1 2 } Pn , and 

/on o — n\d /on 9 — n\d 

\N(V n , mi \ {0}) - { - 4 I < 1 o Pn, 



and 



/on 1 o— n\d /on 1 o— n\ci 

A n R ^ : |iV(l4^\ {0}) - [A Y I < 2 Pn ' ^ 



(2™ + 2" n ) a! (2 n + 2 -nVi 



l^(K, m2 \ {0}) - ^ 2-^-1 < ^ 2~^ Pn 

for all n = r + 1, r + 2, . . . , m x = 1, . . . , (2 2n - m 2 = 1, . . . , (2 2n + l) d , 
Pn = Pn(R) to be defined later. 

First, we establish an upper bound on (1 — T?'(Ar)). The distribution of , 
N(U njmi \{0}) and N {V n ^ mi \{0}) according to the Palm version of the Poisson 

point process is Poisson with mean ^ ^ — '—, and that of N(U nym2 \ {0}) and 



2 

N(V njm2 \ {0}) is Poisson with mean - +2 — — . Thus we can use Lemma [2TT1 
(noting, that p n < 1 for all n > r ) and a simple union bound to get: 



1-P'(A R )< J2 (2 2n -l) d (l-P / (^ 1 ))+ (2 2n + l) d (l-P'(A™' 1 )) < 

n=ro+l n=ro+l 



oo 



2 2 ^ d 8exp( 



2o(n-l)d 
N ,)2(n+l)dg ovri f Pn^ 



n=ro+l 
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Now set p n = 2~ rflf 10" 2 . Then 

1-P'(A R )< J2 2 2 ^ d 8eM~^ 1)d - ^—)< 



n=ro+l 

Rho- 



< ^ 2 2(n+1)fl!+3 exp(-2 (n+1)((i ~t ) 



2 3d 

n=r Q + l 

Observe that the first term in the last series is at most R 2d 2 2d+S exp(—R d ^^-); 
and that for R > Ro(d) large enough (and by that r large enough) the last 
series can be dominated by a geometric series with common ratio smaller 
than 1, and therefore the sum is at most constant times the aformentioned 
first term. Hence there exist constants C, d > 0, such that for every R> Rq 

1 -P'(A R ) < C'exp{-c'R d ). 

Now let us assume that n > tq throughout the following computation. 
Let Q be an arbitrary cube of sidelength 2™ containing 0, U its left and V 
its right side. Then, when conditioned on A R and being a configuration 
point, the following is true: 



N(U) - N{V) 



N(U)+N(V) 



< 4p„ 



We can show this by making an easy observation: there is a cube Q n , mi of 
sidelength 2™ + 2~ n with vertices in G n (y) so that U C U n>mi = U, and a cube 
Qn, m2 of sidelength 2 n — 2~ n with vertices in G n (v) so that U njm2 — U_CU; 

and analogously the same is true for V (note, that U and V are not necessarily 
two halves of the same cube). Note that 

N(U) - N(V) N(U\{0})-N(V\{0}) + 1 



N(U) + N(V) ~ N(U \ {0}) + N(V \ {0}) N(U \ {0}) + N(V \ {0}) ' 

The summands can be bounded using monotonicity of N(.) and the bounds 
provided by A R by its definition : 

N(U\{0})-N(V\{0}) N(U\{0})-N(V\{0}) 
N(U \ {0}) + N(V \ {0}) N(U \ {0}) + N(V \ {0}) 
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< 



^^(1 + P„)-^F^(1-Pn) 1 



2^f^(l-p„) 2 
We make the following remark about p n = p n (R): 

10 -2 2 -(n-ro-l)| < ^ < 10 -2 2 -(n-ro-2)|_ 

Now 



i + Pn / 2" + 2-" y i 

l-p n V2 n -2- n / 



2™ -2- 



1 



A 1 A ' 1 
< ( 1+ 2^) <1 + Crf 2^ 



2 2 ™ - 1 , 

for n positive, and c<j only depends on d (for example q = 2 d ); and 

1 + Pn - 2pr, 



1 - Pn 1 - Pn 

as p n < |, SO 

N(U \ {0}) - N(V \ {0}) 



< 1 + 4p n , 



7V([/ \ {0}) + 7V(V \ {0}) 2 
The second term is bounded by 
1 



< -(4pn + Cd^Z^ + c d ^^A Pn ). 



l 2 2n ~ 



< 2- 



iV(i7 \ {0}) + N(V \ {0}) (2 n -2- n ) d 



1 - p n ) < c d 



2 2n-2 



If i? > i?o(d) for R is suitably large, and therefore p n suitably greater 
than 2 2n ~ 2 for n > r , then symmetry implies that 



N(U) - N(V) 



N(U)+N(V) 



< 4p n 



Let v' be an arbitrary point in M. d , and consider the w'-partition sequence. 
Let u be a point in the smallest cell of 0, and let D n = D n (u,v') denote 
its signed shift in the n'th stage of the AKT(V) along the first axis. Then, 
by the observation that until the n'th stage every displacement takes place 
inside a cube of sidelength 2 n and by the choice of r , we have 



ro 



R 



n=—k 
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if we start the algorithm from the level of cubes of sidelength 2 , where k 
is some non-negative integer. For n > r , D n = (1 — C n {u)) 2 n ~ 1 Mp n ]+N^v] ' 
where Q n is the v' + 2 n Z grid cube containing 0, U n its left and V n its right 
half; and < C n {u) < 1 is the relative distance of u to the bisector of Q n 
orthogonal to the first axis before the n'th stage. Conditioned on Ar and 
being a configuration point we have then: 

n=r +l n=r +l v y v ' n=r + l 

oo oo 

2"4+ 1 j Riio- 2 = 2i?U(r 2 2- 4 2^2~4 < 



=T'o + l n=0 

-1 1 

2"' 



" 4 



since ^ < 2 ro+1 by the choice of r§. Thus we have that, conditioned on the 
event Ar and being a configuration point, for every v' G M d and every w in 
the smallest cell containing in the v ' partition sequence (and in particular, 
for u = 0) the total shift of u along the first axis is at most R. Additionally, 
we get that on Ar the length of the shift in the n'th stage for n > ro(R) 
along the first axis is at most c' n (R) := 2~% +1 RilO~ 2 . For convinience, let 
c' n (R) := 2 n_1 for —k < n < tq, which is the trivial upper bound for \D n \. 
The constructed sequence c' n (R) is absolutely summable, since k e Z + is 
finite almost surely. 

For 1 < i < d one can define the events Ar^ similarly, to bound lengths 
of the shifts along the i'th axis. We use the same p n series, and the bounds 
on the probability of the complement events and lengths of shifts hold here, 
too. 

Now define 

d 

E R :=f)Af ti . (2.1) 

i=l 

This satisfies the conditions of the lemma, and the proof is finished. 

□ 

Fix R > Rq and the corresponding event Er from Lemma 12.21 The 
following lemma tells us that the measure of the cell assigned to by AKT(v) 
run up to stage n tends to 1 with n. This is yet a simple consequence of the 
fact that the number of Poisson points in a cube is concentrated around the 
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volume of the cube. The only extra technicality is coming from the fact that 
we want to prove convergence of the cell volumes uniformly in v, but this can 
be done the same way as in the proof of Lemma 12.21 (using approximating 
cubes that belong to some countable family), hence we omit the details. 

Lemma 2.3. On from Lemma [2~2} the cell C v<n ofO according to AKT(v ) 
satisfies 

lim Vol(C vn ) = lim / f vn (x)dx = 1, 

n— >oo n— >oo Jjjd 

where convergence is uniform in v. 

Next we consider the sidelengths of the cells during the procedure, to 
have a full control on the diameter of C Vjn U {0}. 

Definition 1. Let l(C Vin ) G M. d denote the vector, whose i'th coordinate is 
the i'th sidelength of the cell C v>n . Let |.| denote the maximum norm and 
f{C Vi n) denote maximal volume of the (d — 1) -dimensional faces of C v<n . 

On Er, in each considered cuboid there are at least half as many, and 
at most twice as many configuration points as expected (i.e. the volume of 
the cuboid). So, by the fact that after each AKT(f) step the cells in the 
currently considered cuboids have equal volume, we have 

1 d 

- <\{li{C v , n ) <2, (2.2) 

for all n > r . So for the (d — l)-dimensional volume of the faces of C VtTl we 
have f(C vn ) < —. — — r. 

v v -, n > mini ii(Cv,n) 

Proposition 2.4. Conditioned on Er as in Lemma \2. 21 \l(C Vjn )\ is bounded 
from above by KR uniformly in v, where K > is a constant. 

Proof. It is clear that |/(C„ jn )| < 2 r ° < -j for n < ro, because transforma- 
tions of the cell only happen inside a cube of sidelength 2 r °. 
For n > ro(-R) 

|/(C„, n+1 )| < \l(C v>n )\(l + 4 Pn+1 ), 

because in each stage one sidelength can change in at most one step, and in 
each step of the n'th stage a sidelength can change at most to (1 + 4p n+1 ) 
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times its size on Er. This implies, that 

oo D oo 

\l(C v , n )\ < \l(C v , ro )\ J] (l + 4 Pi ) < -exp( 4 A) < 

i=ro+l i=r +l 

<fexp(| J2 TA Pi )<^ex P (^R)<KR, 

i=ro+l 

for some universal constant K > 0. 

□ 

Proposition 2.5. Conditioned on E R , \l(C VjTl+ i) — l(C Vin )\ < e n (R) forn > 
r (R), where {e n (R)} is an absolutely summable real sequence with sum at 
most cR. 

Proof. From the previous proposition we know that |/(C^ )ri )| < KR, which 

implies that all components of the vector are at most KR. In each stage 

all sidelengths can change at most once. As we are in Er, this change is 

within factors (1 — 4p n ) and (l + 4p n ) in the n'th stage. Then it follows, that 

|/(C„ jn+ i) — l(C v>n )\ < KR4p„. The {p n } series is (absolutely) summable and 

X^ro+i KR^Pn < cR, for some c < constant, which proves the claim. 

_ „ ^ 

Lemma 2.6. There exist c, C > such that for all R > there exists an 
event Er such that 

(i) 1 - P'(Er) < Cexp(-cR d ). 

(ii) Conditioned on Er, the diameter of the cell C v>n of in the AKT(v) 
ran up to stage n is at most c'R, for all v G M. d . 

(tit) On E R we have \\f %n - f v , n +x\\i = Vol(C v , n AC v;n+1 ) < c n (R), for 
n > r Q (R), where the series c n (R) is absolutely summable. Again, the 
constants do not depend on v. 

Proof. We take the same events Er as in Lemma 12.2} therefore (i) is 
satisfied. 

(ii) will be satisfied by the following upper bound on the diameter using 
Lemma [2.21 and Proposition 12.41 

oo 

diam(a, n U {0}) < d\ ^2Di(0,v)\ + d\l(C v , n )\. 

i=—k 
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To verify that (Hi) holds, let us fix n > r (R). Consider the transforma- 
tions of the cell of occurring during the steps of the n'th stage. There are 
d steps, and all can be treated similarly, so we only consider the step along 
the first axis. Let A denote the cell before this step, and B thereafter. We 
introduce an auxiliary cell C: if w G R d is the shift vector of in this step, 
then C := A + w. We bound first Vol(AAC). 

Vol(AAC) < 2\w\f(A) < Pn T+\KR) d - 1 = c'^R), 

where we used the assertion about the maximal volume of a face in Proposi- 
tion [23] and the bound on the length of shifts of provided by the definition 
of En (see the proof of Lemma I2~2l) . Now we bound Vol(CAB) using Propo- 
sition 12.51 

Vol(CAB) < f{A)\l(A) - l(B)\ < (K RY' 1 e n (R) . 

As both bounds considered as series are absolutely summable, it is now easy 
to obtain c n (R) summing the volume differences of the d steps of the consid- 
ered stage and using triangle inequality. 

□ 

3 The allocation rule 

Let c n (R) be as in Lemma [2.61 (Hi). Then we have c n (R) > and the sum 
of the c n (R) over n is finite. Let Er be as in the same lemma. 

Recall that f£ n = f v>n was defined as the indicator function of the cell 
of the origin, where the cell is the result of AKT(f) run up to stage n, and 
u)' G Q' is a configuration containing the origin. 

The next lemma is straightforward from the definitions. 

Lemma 3.1. The map ip : u — > L 1 (M a! ), = jf"^+£ is a weak allocation. 

Lemma 3.2. For P 1 '-almost every u>' G Q' , f v>n is ^-continuous in v, except 
for possibly the union of finitely many hyperplanes. 

Proof. A point of discontinuity v is necessarily such that is on some hyper- 
plane H orthogonal to an axis defined as H = v+{(xx, . . . , Xj-i, i2 e , Xj + i, . . . , Xd)} 
for some j G {1, . . . , d}, i G {1, . . . , 2 n ~ 1 }, and I G {m, m + 1, . . . , n} with 
some m < 0. This m is the biggest nonpositive integer such that inside the 
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cube of in v + 2 n Z d every cube of v + 2 m Z d contains at most 1 configuration 
point. It is easy to check that almost surely there is an m G Z with this 
property that works for every v. 

□ 

Recall from (jl.3p the definition of the averaging function f n — f£ := 
7^3 Jj 2 „] d f v , n dv. The integral in the definition exists by the previous lemma. 

Proposition 3.3. With probability 1, the L 1 limit f of f n exists, it is a 
function with values in [0, 1] of integral 1, and with P'[diam({0}Usupp(/)) > 
r] < Cexp(— cr d ). 

Proof. Fix R, and suppose that u>' G Er. If we prove the proposition for 
R, and show that on E R the diameter of supp/ is bounded by cR with some 
c independent of R, then the bound in the proposition follows. Note that E R 
is a monotone increasing family of events that exhausts a subset of measure 
1 of Q'. Moreover, it is clear from the definition that the function f defined 
for an uj' G Er does not depend on R (so we can define it automatically for 
any element of the increasing union \Jr =1 Er). 
For n < n, u G 2 ft • {0, 1, . . . , 2 n ~ h - l} d } let 

fn,u,u>' fh,u J_ I r 7 

in Jn ■ 9nd / Jv,nUU. 

z Ju+{0,2™] d 

In particular, f n = f™' u . Then 

fn= 2(n-n)d ^ ' 

uG2 fi {0,l,...,2"- fl -l} d 

If u, v! G 2" • {0, 1, . . . , 2 n ~ n - then ffh* = fj' u ' by flE2) (the two inte- 
grated functions are the same). We have \\fn ,u -fl' u \\i = | \2~ nd j u+[0 ^ ]d f v , n dv- 
2 ~ nd I u +[o,2n]d fv,ndv\\i < 2~ M J u+[0)2 a]d \ \fv,n - fv,n\U d v, and similarly for u' . 
By Lemma [2761 (m) this latter is bounded by Y^Ln+i c «(-^)- 
Thus the triangle inequality gives 

oo 

n/r-/r'iii < ii/n ft,u -/?'iii+ii/r-/r1ii+ii/r'-/r'iii <E 2c *w 

(3.2) 
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Using (13.1 p this implies 



Wfn ~ f*'°\\i < £ ^=^ll/r-/: ,0 Hi<2£Q(it!). (3.3) 

«G2 fi {0,l,...,2 n -"-l} d «=n 

On the other hand, for any m,n G Z + , m < n, 

\\f /II ^ \ \ f -f m iO|| i 1 1 fm.O / II 1 1 f /m,0|| i 1 1 fm.O xm,0|| 

UJn- Jm\\l ^ Win- J n lll + ll/n ~ Jm\\l ~ \\Jn ~ J n \\l + \\j n ~ Jrn 111- 

The first term on the right is < 2^°^ m Cj(.R) by (13. 3j) . the second term 
is < Y^i= m c i(R) ^ Lemma [2761 (in). We conclude that (/„) is a Cauchy 
sequence, and so there is a limit / in L 1 . The fact that / takes values in 
[0, 1] follows directly from the same fact to f n . By Lemma 12.31 it is easy to 
see that J f n —>l, and the support of each f n is within radius c'R around 0, 
hence the dominated convergence theorem implies that J f = 1. The above 
hold for every R, hence the proposition follows. 

□ 

The next proposition implies Theorem 11.11 

Proposition 3.4. The family + £ : £ G u} is a fractional allocation 

rule. It satisfies 

P[diam({0} U supp/ w ) > r | G u) < exp(-cr d ). 

Proof. By Proposition 13.31 we have that {f w ~^ + £ : £ G u} satisfies (1') 
and (3) in the definition of a fractional allocation rule. Similarly for the claim 
about the support of / w . Measurability and equivariance are clear from the 
construction. So it only remains to prove (2'). 

Fix R and u G E R . By Lemma EH {f^ + £ : £ G lo} is a weak alloca- 
tion; in particular, (2) holds. Hence for almost every x G M d , Yl^ifn" + 

00*) = E^ nd U [0 ,2^fv^Hdv)(x) = 2- nd J mn]d (Z^f^+h))dv = 

2~ nd L 2 „ ]d ldv = 1, snowing that + £ : £ G w} satisfies (2'). Finally, 

the fact that ^£ gw + £ is identically 1 almost everywhere implies this 
for E €e . limn /r € + e- 

□ 

Lemma 3.5. Let ,x G lR d . T/ien P-almost surely there is an e and £i, • • • , G 
u; stzc/i #icrf t/ie e-neighborhood B x (e) of x is contained in nj =1 supp(/ w_ ^+£ i ) ; 
and B^e) n suppC/""* + = for every (Gw \ {6, . . . ,&}. 
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Proof. Let e' > be arbitrary. For almost every uj the sum 

P[diam(supp(r^)) > \x - £\ + e'] 

few 

is finite, as a consequence of Lemma 12.61 Hence almost surely only finitely 
many of these events happens, by the Borel-Cantelli lemma. 

Thus B x (e') is intersected by finitely many of the supp(/ w ~'») (£ G uj), and 
for small enough e > and some £i, . . . , G uj, B x (e) C nf =1 supp(/ w ~^ +6)- 

□ 

A direct consequence is the following: 

Lemma 3.6. For P-almost every uj one can partition M. d to countably many 
path- connected measurable sets of the form S = nf =1 supp(/ aJ ~^ + £*) (with 
some £i, . . . ,£ k G wj. 

The previous lemma enables us to define an allocation rule from our 
fractional allocation rule, in such a way that the cell allocated to a center 
£ G u; is contained in supp(/ w ~^+£). Namely, for each set S as in Lemma [3T6| 
partition S into measurable pieces Si, . . . , Sk such that X(Si) = J s f^~^ + 
£(x)dx. Do it in a way such that S determines the pieces Si in some previously 
fixed (deterministic) way, and such that the pieces change continuously with 
S (in terms of Hausdorff distance between sets, say). A method to do so was 
suggested to us by Yuval Peres, replacing the original, less elegant proof for 
the following lemma: 

Lemma 3.7. Let {f u ~^ + £ : £ G uj} be a fractional allocation rule that 
satisfies Lemma WTSi Then there is an allocation rule : £ G oj} such 

that for every £ G uj we have supp(^(£)) C supp(/ w ^). 

Proof. For each set of the form S = nf =1 supp(/ w ^ ! + £j) as in Lemma [3761 
let q := J s (/ £J '~^+£i)(x)dx, i = 1, . . . , k. So Yli=i c « ^ s ^ ne Lebesgue measure 
of S. Apply a version of the Gale and Shapley algorithm (see [1] for a more 
detailed description) within S as follows. Let each £j start growing a ball 
around itself at linear speed, simultaneously. At time t let each £j capture all 
points of S that its growing ball meets, as long as no other ball has captured 
the point at some earlier time. However, when the total measure of points 
in S that & has captured by time t is equal to q, stop growing its ball, that 
is, if Pi(t) is the set of points in S that have been captured by £j by time t, 
let Pi(t') be equal to Pi(t) as long as t' > t. It is intuitively clear that by 
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some time t the Lebesgue measure of Pi(t) is equal to q for each i. Define 
this Pi{t) to be the part of the cell of £ within S. Doing this for each of 
the countably many S, we get the cell of & assigned to it by our allocation. 
Measurability and invariance follow from our method and the assumptions; 
see [1] for more details and rigorous arguments. 

□ 

We finish by noting that Proposition 13.41 and Lemma 13.71 imply Theorem 

If we wanted the cells assigned by our allocation rule to be connected and 
to contain the corresponding centre, this could be done by growing "tendrils" 
that connect the pieces of each cell C of the original allocation, and the center 
for C . With care about preserving measurability and equivariance this can 
be done, but we omit the details here. 

Remark 3.8. An intuitive interpretation of the AKT algorithm comes from 
thinking of the configuration points as gas particles. Then the procedure 
systematically equates pressures between neighboring cubes in the dyadic 
subdivision. For the sake of analysis it was easier to allow cells around par- 
ticles to expand in the directions of the axes and by, but a more "canonical" 
version would be obtained without this artifact. That is, put one particle 
of gas in each center of the point process, and start by growing cells around 
them (small balls at the beginning), whose pressures would be proportional 
to the volume of the cell and at normal direction to the surface. When two 
cells meet, the pressure differences between them would tend to equate, and 
in the limit all the cells would have the same volumes. The tail behavior of 
their diameters should be as good as that of our fractional allocation rule 
(and perhaps better if we take constant factors into account). On an even 
more speculative note, we mention that the above procedure looks like a 
modification of the stable allocation rule: do not fix centers, and let the 
growing cells "push" each other while occupying yet unoccupied territories. 

Question 3.9. Can one make the above heuristics precise, to obtain a canon- 
ical allocation rule of optimal tail? 
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